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LONG-TERM  GOALS 

This  work  seeks  to  understand  the  effects  of  topography  on  the  ocean  general  and  regional  circulation 
with  a  focus  on  the  wide  range  of  scales  of  interactions.  The  small-scale  details  of  the  topography  and 
the  waves,  eddies,  drag,  and  turbulence  it  generates  (at  spatial  scales  ranging  from  meters  to 
mesoscale)  interact  in  the  boundary  layers  to  influence  the  larger-scale  flow.  We  will  study  these 
issues  through  ocean  model  simulations,  adjoint  sensitivity  experiments,  and  state  estimation  using 
measurements  in  the  region  surrounding  an  island  in  the  westward-flowing  limb  of  the  subtropical 
gyre. 

OBJECTIVES 

The  objectives  of  this  program  include  addressing  number  of  questions  on  ocean  circulation  that  arise 
near  islands  and  abrupt  topography.  Questions  to  study  include: 

•  What  is  the  steady-state  flow  around  the  feature,  including  isopycnal  depths,  currents,  and 
water  mass  characteristics?  How  do  these  vary  with  changes  in  the  surrounding  large-scale 
flow  and  atmospheric  forcing? 

•  What  is  the  response  of  the  local  flow  to  time-dependent  large-scale  flow?  i.e.,  how  does 
the  feature  scatter  incoming  barotropic  and  baroclinic  Rossby  waves  (as  well  as  near- 
inertial  waves  and  tides),  into  outgoing  waves,  mean  flows,  and  topographically-trapped 
and  damped  waves? 

•  How  much  can  be  learned  about  the  regional  circulation  or  the  larger  scale  currents  from 
local  observations?  Can  they  be  used  for  linear  or  nonlinear  inverse  estimates?  i.e.,  can  the 
regional  or  gyre  circulation  be  monitored  from  pressure  gauges,  temperature  sensors, 
current  meters,  or  other  measurements  near  the  feature? 

•  The  influence  of  the  boundaries  on  the  flow  in  the  nested  inner  model  can  also  be  turned 
around  to  study  how  the  feature  influences  the  larger-scale  boundary  flow.  Example 
questions  are:  what  fonn  drag  is  exerted  on  the  larger  scale  flow?  How  far  eastward,  or  “up- 
current”  does  the  influence  of  the  feature  extend?  How  turbulent  are  the  wakes?  Are  eddies 
shed?  If  so,  do  shed  eddies  create  counter-currents  in  the  wake  of  the  feature? 
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APPROACH 


We  will  explore  and  answer  these  questions  through  numerical  model-based  data  analysis,  including 
state  estimation  in  nested  domains.  We  will  use  assimilated  global  HYCOM  and  the  MIT  general 
circulation  model  (MITgcm)  in  nested  domains  starting  from  basin  scale,  connected  to  high-resolution 
(sub  km)  inner  domains  modeled  with  Delft3D  or  ROMS  within  about  a  degree  of  the  island  which 
will  be  configured  and  run  by  other  investigators.  Our  results  will  be  used  in  collaboration  with  other 
investigators  to  plan  the  experiments  and  to  help  to  understand  the  observations.  We  will  use  the  ocean 
models  to  provide  a  dynamically-consistent  re-analysis  of  the  observations  in  the  region,  both  from  the 
experiment  and  from  other  sources  (Argo,  Altimetry,  SST,  TAO,  etc.). 

WORK  COMPLETED 

We  have  been  investigating  model  representations  of  the  interactions  of  the  NEC  with  the  archipelagos 
and  ridges  in  the  western  Tropical  Pacific  using  analyses  from  both  global  HYCOM/NCODA  and 
MITgcm  state  estimates.  This  is  an  area  where  topography  varies  much  faster  than  allowed  by  Quasi- 
Geostrophic  theory,  and  it  is  not  clear  how  well  the  large-scale  models  adequately  represent  the 
physics  of  the  topographic  interactions.  In  addition  to  the  HYCOM  output,  we  are  using  a  complete 
Tropical  Pacific  model  (1/3 “horizontal  resolution,  30°S  to  30°  N)  as  well  as  a  regional  1/6°  model 
originally  developed  for  the  ONR  project  "Origins  of  the  Kuroshio  and  Mindanao  Currents"  (OKMC). 
Both  models  include  state  estimation.  The  basin  model  resolves  wind-forced  currents  propagating 
westward  approximately  along  Rossby  wave  characteristics  as  can  be  seen  by  calculations  of  the 
sensitivity  of  western  boundary  transport  to  wind  stress  [Zhang,  et  al.  2012].  A  comparison  of  the  SSH 
estimates  from  2010-2011  from  AVISO,  a  free  forward  simulation  of  the  OKMC  regional  model  using 
HYCOM/NCODA  initial  and  boundary  conditions,  and  Tropical  Pacific  and  OKMC  regional  state 
estimates  using  existing  data  (such  as  Argo,  SST,  and  SSH)  is  shown  in  Figure  l(a,b,c)  for  three  1°  by 
1°  boxes  just  to  the  east  of  Palau,  Yap,  and  Guam,  respectively.  Although  the  AVISO  map  is  only  an 
estimate,  it  is  taken  as  truth  for  these  comparisons.  The  state  estimates  used  along-track  SSH  from 
three  satellites  (Jason- 1,  Jason-2  and  Envisat)  from  the  RADS  database,  so  they  are  not  required  to 
match  mapped  AVISO  SSH  exactly.  The  free  model  simulation  reproduces  some  of  the  SSH 
variability,  but  often  has  quite  large  errors.  The  OKMC  estimate  is  best  for  Palau,  which  is  well  inside 
its  domain  of  interest,  but  the  Tropical  Pacific  estimate  is  generally  better  for  Yap  and  Guam,  which 
are  near  the  eastern  boundary  of  the  OKMC  regional  model,  and  do  not  change  much  from  the 
HYCOM  solution  due  to  the  influence  of  the  HYCOM  boundary  conditions.  The  AVISO  estimates 
show  similarities  between  Palau  and  Yap  (Figure  Id),  which  is  expected  since  Yap  is  just  to  the  east  of 
Palau,  and  westward -propagating  features  should  show  up  similarly  at  both  locations.  In  this 
comparison,  the  phase  lag  seems  to  be  about  1  month,  but  there  are  significant  differences,  due  to  local 
wind  stress,  observation  error,  or  alternative  pathways. 

The  state  estimates  will  be  used  for  diagnosis  of  the  physical  processes  active  around  the  features  after 
comparisons  with  independent  observations,  testing  the  model  hindcasts  and  forecasts.  For  example, 
the  model  performance  in  the  highly-stratified  positive-wind-curl  region  is  quite  sensitive  to  the 
assumed  vertical  diffusivity  coefficient.  We  are  taking  guidance  from  estimates  of  the  diffusivity  from 
Argo  (Whalen,  Talley,  and  Mackinnon,  2012)  and  refining  the  background  diffusivity  estimates  for 
this  region.  The  state  estimates  are  being  re-run  with  other  adjusted  parameters,  such  as  viscosity,  and 
with  varying  assimilation  windows  to  see  how  long  an  assimilation  window  can  be  used  and  how  well 
the  estimates  perform  in  forecasts. 
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In  addition,  MITgcm  simulations  of  an  idealized,  flat-bottom,  rectangular  ocean  basin  containing  an 
idealized  island  are  being  used  to  study  the  sensitivity  of  transport  to  wind  stress  in  the  vicinity  of  the 
island  and  through  the  section  connecting  the  island  to  the  eastern  boundary,  generalizing  the  island 
rule.  This  is  being  compared  to  simpler,  semi-analytic  island  rule  streamfunction  calculations, 
although  the  idealized  calculations  are  surprisingly  complicated  and  hard  to  get  right  (Kessler,  personal 
communication). 

RESULTS 

The  waveguides  present  in  the  basin  (through  beta  and  topography)  determine  the  pathways  of  control 
of  the  transport,  and  are  hypothesized  to  correspond  to  the  integration  contour  for  the  island  rule 
calculation,  although  sharp  dependence  on  a  single  latitude  is  not  observed. 

IMPACT/APPLICATIONS 

The  boundary  current  transport  sensitivity  to  local  and  basin-wide  wind  stress  is  quantified  for  time- 
dependent  cases  with  arbitrary  topography  and  existing  background  flows. 
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SSH  comparison:  Palau 


SSH  comparison:  Yap 


Figure  1  a,b,c:  Comparison  of  AVISO  mapped  SSH  (blue  line)  with  a  free  1/6  degree  OKMC 
regional  simulation  (red  line),  a  Tropical  Pacific  state  estimate  made  over  2  month  intervals  (green 
line),  and  a  state  estimate  for  the  OKMC  regional  model  using  one  month  assimilation  windows 
(cyan).  The  values  are  averages  in  1°  by  1°  boxes  just  to  the  east  of  Palau,  Yap,  and  Guam, 
respectively.  Id  is  a  comparison  of  the  AVISO  SSH  for  Yap  (blue)  and  Palau  (green),  repeated  from 

figures  la  and  lb. 
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